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Abstract 

Strong consistency of the maximnm likelihood estimator (MLE) for parametric Gibbs point 
process models is established. The setting is very general. It includes pairwise pair potentials, 
finite and infinite multibody interactions and geometrical interactions, where the range can be 
finite or infinite. The Gibbs interaction may depend linearly or non-linearly on the parameters, 
a particular case being hardcore parameters and interaction range parameters. As important 
examples, we deduce the consistency of the MLE for all parameters of the Strauss model, the 
hardcore Strauss model, the Lennard-Jones model and the area-interaction model. 
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1 Introduction 


They are now 
image 


Gibbs point processes are popular and widely used models in spatial statistics to describe the 
repartition of points or geometrical structures in space. They initially arose from statistical physics 
where they are models for interacting continuum particles, see for instance 2 
used in as different domains as astronomy, biology, computer science, ecology, forestry 
analysis and materials science. The main reason is that Gibbs point processes provide a clear 
interpretation of the interactions between the points, such as attraction or repulsion depending 
on their relative position. We refer to i, [II, 113 and for classical text books on spatial 
statistics and stochastic geometry, including examples and applications of Gibbs point processes. 

Assuming a parametric form of the Gibbs interaction, the natural method to estimate the 
parameters is likelihood inference. A practical issue however is that the likelihood depends on an 
intractable normalizing constant, called the partition function in the statistical physics literature, 
that has to be approximated. Some sparse data approximations were first proposed in the 80’s, 
e.g. in [^ . before simulation-based methods have been developed 111. A comparative simulation 
study carried out in demonstrates that Monte Carlo approximation of the partition function 
provides the best results in practice. To avoid the latter approximation, other estimation methods 
have been introduced, including pseudolikelihood and moments based methods, see the books 
cited above. With modern computers, the Monte-Carlo approximation of the partition function is 
no longer an important issue and maximum likelihood estimation for spatial data is feasible and 
widely used in practice. 

From a theoretical point of view, very few is known about the asymptotic properties of the 
maximum likelihood estimator (MLE) for Gibbs point processes. The asymptotic here means that 
the window containing the point pattern increases to the whole space It is commonly believed 
that the MLE is consistent and more efficient, at least asymptotically, than the other estimation 
methods. This conviction has been supported by several simulation studies, ever since (2^, see 
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also The present work is concerned with the basic question of consistency. The latter is con¬ 
jectured to hold in a very general setting but no proof were so far available in the continuous case 
of Gibbs point processes on This is in contrast with the discrete case of Gibbs interactions on 
a lattice, where consistency is established for most standard parametric models, regardless of the 
occurrence of phase transition (when the Gibbs measure is not unique), see [l^ and 121. In fact, 
the continuous case is more challenging in that most parametric Gibbs models involve irregular 
parameters as hardcore parameters (controlling the support of the measure) or interaction range 
parameters. These specificities result in a discontinuous likelihood contrast function, even asymp¬ 
totically, and some further technical difficulties (for instance the true unknown Gibbs measure is 
not absolutely continuous with respect to the Gibbs measure associated to an estimation of the 
hardcore parameter). In an unpublished manuscript [l8l |. S. Mase addressed the consistency of the 
MLE for superstable and regular pairwise interactions (a formal definition will be given later). His 
main tool was the variational principle for Gibbs processes, following the initial idea developed in 
the discrete case in 1^. He restricted his study to pairwise interactions that are linear in their 


parameters, which yields a convex contrast function. His result do not imply any restriction on the 
parameter space, thus including the possibility of phase transition, but the parametric interactions 
considered in [l8| remain nonetheless rather restrictive and do not include hardcore or interaction 
range parameters. 

We prove in this paper that the MLE is strongly consistent for a wide class of stationary Gibbs 
interactions, without any major restriction on the parameter space. Our assumptions include finite 
and infinite-body interactions with finite or infinite range, and we do not assume any continuity 
with respect to the parameters. Our result covers in particular the consistency of the MLE of all 
parameters of the Strauss model (including the range of interaction), the hardcore Strauss model 
(including the hardcore parameter), the Lennard-Jones model (including the exponent parameters), 
and the area-interaction process (including the radius of balls). An important ingredient of the 
proof is the variational principle, as in |16l | and 181, which guarantees the identifiability of the 


parameters. Our original contribution is the formulation of a minimum contrast result in presence 
of discontinuities and, in order to apply it, new controls of physical quantities as the pressure and 
the mean energy with respect to the parameters. 

Beyond consistency, the next natural question concerns the asymptotic distribution of the MLE. 
This problem is even more arduous and we do not address it in the present paper. Nonetheless, let 
us briefly mention the state of the art on this question. In 1992, S. Mase [l7| proved that the MLE 
of regular parameters in a certain class of Gibbs interactions is uniform locally asymptotic normal. 
However his result, relying on strong cluster estimates, is established for Gibbs measure generating 
very sparse point patteri^ which implies restrictive conditions on the parameter space. In the 


same period, J. Jensen |14| proved the asymptotic normality of the MLE under the Dobrushin 
uniqueness region, where the Gibbs process satisfies mixing properties. Here again, as noticed in 
[HI , this assumption implies strong restrictions on the parameter space. Without these conditions, 
phase transition may occur and some long-range dependence phenomena can appear. The MLE 
might then exhibit a non standard asymptotic behavior, in the sense that the rate of convergence 
might differ from the standard square root of the size of the window and the limiting law might be 
non-gaussian. To our knowledge, the only well understood example is the estimation of the inverse 
temperature in the Ising model on a square lattice studied in (2^. For this example, the asymptotic 
law of the MLE is always Gaussian and the rate of convergence is standard except at the critical 
temperature where it is faster. Generalizing this result to other Gibbs models, especially in the 
continuous case, is hazardous. The main reason is that the occurence of phase transition is in 
general not well understood, and this phenomenon can be of very different nature depending on 
the model. Gharacterizing the asymptotic distribution of the MLE in a general setting of Gibbs 
models still remain a challenging open question. 

The remainder of the paper is organized as follows. Parametric Gibbs point processes and the 
MLE procedure are described in Section [2l Section |3] contains our main result, namely minimal 
conditions on the Gibbs interaction to ensure strong consistency of the MLE. While the latter 
result is established in a very general setting, we present in Section |4] standard families of models 
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where our main result applies. Specifically, we deal with finite range pair potentials with or 
without hardcore (including the Strauss model and the hardcore Strauss model), infinite range 
pair potentials (including the Lennard-Jones model), and infinite-body interactions (including the 
area-interaction model). Section [S] contains the proofs of our results. 

2 Gibbs point processes and the MLE 

2.1 State space, reference measure and notation 

We consider the continuous space of arbitrary dimension d > 1. The Lebesgue measure on 
is denoted by and the symbol A will always refer to a bounded Borel subset of For a; S K'^, 

I a: I denotes the Euclidean norm of x while for A C |A| := A‘^(A). 

A configuration is a subset w of which is locally finite, meaning that w n A has finite 
cardinality Na{uj) := #{uj fl A) for every bounded Borel set A. The space SJ of all configurations 
is equipped with the cr-algebra J- generated by the counting variables A/a- 

For convenience, we will often write wa in place of wfl A and wa'^ for wfl A'^ = uj\u}\. Similarly 
for every w and every a; in w we abbreviate w U {a;} to w U a; and a;\{a;} to U!\x. 

As usual, we take the reference measure on (fJ, iF) to be the distribution tt of the Poisson point 
process with intensity measure A'^ on Recall that tt is the unique probability measure on (12, if) 
such that the following hold for all subsets A: (i) A^a is Poisson distributed with parameter A‘^(A), 
and (ii) given N\ = n, the n points in A are independent with uniform distribution on A. The 
Poisson point process restricted to A will be denoted by tta. 

Translation by a vector u S is denoted by Tu, either acting on or on 12. A probability P 
on 12 is said stationary if P = P o for any u in R'^. In this paper we consider only stationary 
probability measures P with finite intensity Ep{N[Q ijd), where Ep stands for the expectation with 
respect to P. We denote by P the space of such probability measures. 

We denote by Ag, P„ and A„ the following sets 

Ao = [0,1)"^, Xn = and = [J Tfc(Ao) = [-n, 

keXn 

In the following, some infinite range interaction processes will be considered. To ensure their 
existence, we must restrict the set of configurations to the so-called set of tempered configurations 
as in [ 2 ^. The definition of the latter may depend on the type of interactions at hand. In case of 
superstable pairwise interactions, it is simply defined by 

12t = {w S 12; 3t > 0,Vn > 1, ^ -^Ti(Ao)(^) — i(2n)‘^}. 

iein 

Unless specified otherwise, this is the definition we consider in the sequel. From the ergodic 
theorem, any second order stationary measure on 12 is supported on 12-^) so the restriction to 12 t 
is a mild assumption in practice. 

Several other notation are introduced throughout the next sections. For convenience to the 
reader, we summarize the most important of them below. 
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ri, rix 

Wa, Wa<= 

5, e 

/, 0, /c 

HUoj) 

rS,9 / \ 

/a (w) 

rySfi 

gsfi 

G 

LU* 

H\P),I{P) 

p{6,9) 


space of configurations, tempered configurations respectively 

configuration w inside A, outside A respectively 

parameters of the interaction where S is the hardcore parameter 

parameters space: 5 € I = [i5niin, ^max], 0 € Q and /C C 0 is compact 

space of configurations with locally finite energy 

energy (or Hamiltonian) of a; in A 

conditional density of w inside A, given uja<=, see 0 

partition function associated to the free boundary condition, see (HU 

set of stationary Gibbs measures associated to w !->■ f^iuj) 

G = 'JseipGeG^'^ 

configuration associated to the true unknown parameters and 0* 
likelihood contrast function for the observation of to* on A„, see ® 
mean energy and specific entropy oi P G G, see (ESI and (HU) 
pressure associated to P^’^ G G^'^, see m 


2.2 Gibbs point processes 

From a general point of view, a family of interaction energies is a collection 'H = (i?A), indexed by 
bounded Borel sets A, of measurable functions from GIt to K U {+oo} such that for every A C A', 
there exists a measurable function (/3 a,A' from GIt to K.U {+oo} and for every ui G GIt 

Ha'{lo) = Ha{uj) + ( 1 ) 

This decomposition is equivalent to (6.11) and (6.12) in p. 92]. In physical terms, Ha{uj) = 
PIa(u}\ U cua<=) represents the potential energy of the configuration uja inside A given the config¬ 
uration ccao outside A. In words, o is just a compatibility property stating that the difference 
between the energy on A' and the energy on a subset A only depends on the exterior configuration 

UJAc. 

As explained in introduction, we aim at including a large class of interactions in our study, 
to cover for instance infinite-body interactions like the area-interaction process considered in Sec¬ 
tion 14.41 Nevertheless, most of standard parametric Gibbs models are pairwise interaction point 
processes. They are introduced below and we will come back to this important class of Gibbs 
models in Section 01 Pairwise potential interactions take the particular form 

Ha{uj) = z Na{uj) + ^ (j){x-y), (2) 

where z > 0 is the intensity parameter and (p is the pair potential, a function from to RU{-|-oo} 
which is symmetric, i.e. </>(—x) = (f){x) for all x in In connection with ([ij, if Ha satisfies ®, 
then (/ja.aKwaO = z iVA'(wA=) + Y.{x,v}(iu,A<^,{x,y}nK'^<h - y)- 

Let us present two well-known examples. We will use them through the paper to illustrate our 
notation and assumptions. 

Example 1: The Strauss pair potential, defined for some possible hardcore parameter 5 > 0, some 
interacting parameter /3 > 0, and some range of interaction R > 0 corresponds in ([2]) to 

{ oo if jxj < S, 

p if <5 < jxj < R, (3) 

0 if jxj > R, 

if i? > ^, while (p{x) = ooII| 2,|<5 if i? < ^. As we will see later, the theoretical properties of this 
model strongly differ whether <5 = 0, which corresponds to the standard Strauss model, or 5 > 0, 
which is the hardcore Strauss model. 
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Example 2'. The general Lennard-Jones (n, m) pair potential is defined for some d < m < n, and 
some A > 0, i? G M by (m with 

0(x) = A|a:r” - Blccr™, x G (4) 

The standard Lennard-Jones model in dimension d = 2 and d = 3 corresponds to n = 12 and m = 6. 


For a family of energies H — (Ua), we denote by floo the space of configurations which have a 
locally finite energy, i.e. w G floo if and only if, for any bounded Borel set A, H[^{uj) is finite. 

The Gibbs measures P associated to H are defined through their local conditional specification, 
as described below. For every A and every w in floo H Ht, the conditional density /a of P with 
respect to tta is defined by 


/a(w) 


■^a(wa‘=) 


(5) 


where Za{uja<^) is the normalization constant, or partition function, given by 


ZaM = J 


Some regularity assumptions on H are required to ensure that 0 < ZA(a;Ac) < +oo, implying that 
the local density is well-defined. They will be part of our general hypothesis in Section [3] and are 
fulfilled for all our examples in Section 01 

We are now in position to define the Gibbs measures associated to H, see for instance 0- 

Definition 1 A probability measure P on is a Gibbs measure for the family of energies 'H if 
P(f2oo n fix) = 1 and, for every bounded borel set A, for any measurable and integrable function g 
from n to R, 


g{w)P{du}) = J J g{u}ALIu!A‘=)fA{u}AGu;A<^)TTA{dujA)Pidw). 


( 6 ) 


Eguivalently, for P-almost every ui the conditional law of P given wa^ is absolutely continuous with 
respect to tta with the density /a defined in ®. 

The equations (jb]) are called the Dobrushin-Lanford-Ruelle (DLR) equations. Conditions on 
"H are mandatory to ensure the existence of a measure P satisfying ®. For general interactions, 
we will later assume that both (0]) is well-defined and P exists. On the other hand, the unicity 
of P does not necessarily hold, leading to phase transition. As explained in introduction, our 
consistency results are not affected by this phenomenon. 


2.3 Parametric Gibbs models and the MLE procedure 

We consider a parametric Gibbs model that depends on a hardcore parameter J G /, where 
I = [i^min, ^max] with 0 < Jmin < <5max < 00 , and On a parameter 0 G 0 with 0 C and 0^0. 
We stress the dependence on these parameters by adding some superscripts to our notation, e.g. the 
partition function Za{uja'=) becomes Z^’^(wa=) and the interaction energy Ha{uj) becomes H^{u}). 
In particular, the reason why we write Elf^{u}) and not El^pf {uj) (and similarly for other quantities) 
is due to the fact that under our assumptions, the interaction energy will not depend on 5, as 
explained below. 

Specifically we assume that the conditional density writes for any S € I and any 0 G 0 


fS,e 

JA 


M = 


-H 




(w), 


G n 


T, 


(7) 


where 

w G nfo \x — y\> S, 
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while for any 9 G Q and any uj G fir, H^(w) < oo. This specific parametric form clearly indicates 
that the hardcore parameter 5 only rules the support of the measure, and has no effect on the 
interaction energy H^. This is what happens for most models in spatial statistics, see Section H) A 
counter-example is the Diggle-Gratton model Section 6], where the interaction energy depends 
also on 6. This situation is not covered by our study. 

Let be the set of stationary Gibbs measures defined by the conditional density © and 
with finite intensity. A minimal condition is to assume that for any value of the parameters, 0 is 
well-defined and that the latter set is not empty. 

[Existence]: For any S G I, any 9 G 0, any bounded set A in and any cj G fl fir, 
Z^^(uja<=) < oo and the set is not empty. 


Example 1 (continued): For the hardcore Strauss model defined in ([3]) with (5 > 0, the hardcore 
parameter is <5 and the parameter 9 corresponds to 9 = {z,R,p). The assumption [Existence] 
for this model holds if we take I = R+ and 0 = ]R_|_ x K_|_ x R. On the other hand, for the 
standard Strauss model corresponding to the case <5 = 0 in ([3]), there is no hardcore parameter, 
meaning that / = {0} and the only parameter of the model is 9 = {z,R,f5). The existence of 
this model, i.e. [Existence], holds iff 9 belongs to 0 = K+ x K_|_ x K_|.. Note that B needs to 
be nonnegative when there is no hardcore, contrary to the hardcore case. We refer to [26j| for the 
existence results involving stable and superstable pairwise potentials, a class of interactions which 
includes the Strauss model. 


Example 2 (continued): For the Lennard-Jones model defined in ([3]), there is no hardcore param¬ 
eter (J = 0) and the parameter 9 corresponds to {z,A,B,n,m) G 0 C where 0 i s just defined 
via the constraints z > 0, A > 0 and d < m < n. The existence is also proved in 


26|. 


Let 6* G I and 9* in 0. We denote by w* a realization of a Gibbs measure belonging to . 
The parameters and 9* represent the true unknown parameters that we want to estimate. The 
MLE of {6*, 9*) from the observation of w* on A„ is defined as follows. 

Definition 2 Let 1C be a compact subset of 0 such that 9* G 1C. The MLE of S* and 9* from the 
observation of oj* in A„ is defined by 

{L Jn) = argmax ). (8) 

{SM)eIxK 

In Lemma [2] below, we give sufficient conditions which ensure the existence of the argmax in 
([8]). Let us note that we consider here the MLE with free boundary condition, meaning that 
the configuration outside A„ is not involved. This is the most natural setting given that we 
observe w* only on A„. A MLE procedure that depends on the outside configuration could have 
been considered as well and similar theoretical results would have been proved. To implement this 
alternative procedure in practice yet, the interaction has to be finite range so that the outside 
configuration reduces to boundary effects that can be handled by minus sampling. For infinite 
range potentials, as considered in Section 031 this method is not feasible and the MLE with free 
boundary condition makes more sense. 

We will use in the following the equivalent definition of the MLE 

{5n,0n)= argmin 
{S,e)GlxK 


where is the contrast function 




HZ'/J 

|A„| 


Hf 


KJ 


|A„ 




( 9 ) 
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with 


and 


Sn{uJA„) = min \x - y\ 

(10) 

{x,y}euiAn 

= J e-^"^<-‘^^HnsJujA)7CAiduJA)- 

(11) 


In statistical mechanics, the first term in (jH]) is called the finite volume pressure while the second 
term corresponds to the specific energy. 

In the following lemma we give an explicit expression for S„. 


Lemma 1 Under [Existence], the MLE of the hardcore parameter S is 6^ = 6n if € I and 
5n = <^max otherwise, where Sn is given by (fTUl) . consequently 


On = argmin ). 

6»e/c 


( 12 ) 


Proof. If 5 > Sn, then = oo and so Sn < Sn- For any S < Sn and any 0 G 0, 

= |A„|“^(ln(Z^’^) + Note that if 5 < S' , then C and from (ITTI) the 

function S H> is decreasing. Therefore, for any 0 G 0 and for any 5 <5 !->■ is 

decreasing, proving that Sn Sn- Hence Sn = Sn if Sn G I- The other statements of Lemma [T] are 
straightforward. ■ 


This lemma provides a useful result for practical purposes, since it states that the MLE for the 
hardcore parameter is just the minimal distance between pairs of points observed in the pattern, as 
it is usually implemented. Let us remark that a crucial assumption is that the space of parameters 
in [Existence] is the cartesian product 1x0, i.e. the parametric model is well-defined for any 
value of 0 G 0, whatever S is. This is the case for the examples considered in Section H) 

As explained in introduction, the computation of 0„ given by (TT^ is more difficult because there 
is in general no closed form expression for . In practice the optimization (1121) is conducted 
from an approximated contrast function, based on Markov chain Monte Carlo methods where the 
error of approximation can be controlled, see for instance [ill ]. Even to ensure the existence of 0„, 
we need some extra assumptions, as detailed in the following lemma. 

[Argmax]: For any w G and any n > 1, the function 0 i— ^ H^^(uja„) is lower semicontinuous 
over /C, i.e. for any 9 & 1C liminfg/H^g and there exists an upper semicon¬ 
tinuous version 0 >->• H^^{oja„) such that for every 0 G /C, = ^A„ 7 ’'a„- almost surely. 

Example 1 (continued): For all x ^ 0, the function 0 i— >■ 4>{x) in ([3|) is lower semicontinuous which 
ensures that the associated energy 0 i—>■ Hj( for the Strauss model is lower semicontinuous as well. 
Changing the value of (j) at the discontinuity points, it is easy to obtain an upper semicontinuous 
version of the energy function which is 7rA„ -almost surely equal to Therefore [Argm 2 Lx] 

holds for the Strauss model. 

Example 2 (continued): [ArgmaLx] also holds for the Lennard-Jones model since the function 
0 !->■ (f){x) in dH) is continuous for any x 

Lemma 2 Under [Existence] and [Argm2Lx], the MLE of 6 exists, ie. there exists at least one 
On ClC such that = mSngfzfc - 

Proof. It is sufficient to show that 0 >->■ is lower semicontinuous. As a direct 

consequence of [Argmax], 0 i—^ E[f^^{u}\^) is lower semicontinuous. For the finite volume pressure. 
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this comes from [Argmarx] and the Fatou’s lemma since 



3 Consistency of the MLE for general Gibbs interactions 

We detail in this section the minimal assumptions on the parametric family of Gibbs measure that 
imply the strong consistency of the MLE. All of them are fulfilled by the examples of Section |4l 
The first ones, [Existence] and [Argmax] introduced in the previous section, imply the exis¬ 
tence of a Gibbs measure for any 6 G I and 6 G 0, and the existence of the MLE. 

We assume the following mild assumptions on the family of energies. They gather the stability 
of the energy function, a standard stationary decomposition in terms of mean energy per unit 
volume and they deal with boundary effects. We put Q = ■ 

[Stability]: For any compact set /C C 0, there exists a constant k > 0 such that for any A, any 
9 G K. and any lu G ^It 


HiiujA) > -kNa{u}). 


(13) 


[MeanEnergy] : There exist measurable functions Hq and dH^ from fir to K such that for all 
n > 1 and all 0 G 0 



(14) 


where for all P G G and all 0 G 0 the mean energy of P defined by 


H\P)=Ep{H^^) 


(15) 


is finite and for any compact set /C C 0, for P-almost every lo G Glp, 



[Boundary]: For all P gG, for any compact set /C C 0 and for P-almost every uj G Vlp 



To prove the consistency of the MLE, we need to control the asymptotic behaviour of the fi¬ 
nite volume pressure and the specific energy in ([9|). The following assumption deals with the latter. 


[Regularity]: Eor any P in G, for any compact set /C C 0, there exists a function g from R+ to 
K+ with limr_).o gir) = 0 such that for any 6 G K and any r > 0 



( 16 ) 
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Moreover there exists a bounded subset Ag such that for any rj > 0 and any 6q G 1C, there exists a 
finite subset Af{9o) C B{9o, r])r\Q and r{9o) > 0 such that r{9o) < ij and for any A D Ag, 


max inf inf 


( HHuja) - hUuja) 

y Aa(wa) 


> -9{r{do))- 


(17) 


Assumption (1161) implies in particular that the mean energy H^{P) is continuous with respect 
to 9. Let us note that this does not imply that the specific energy or 0 i—>■ Hq is continuous with 
respect to 9. As a counterexample, the Strauss model in Example 1, or more generally any model 
from Sections Q] and EM has a discontinuous specific energy but the mean energy is continuous, 
which is verified in the proof of Theorem [2j Assumption (ITTI) is related to the regularity of the 
coefficient n in [Stability] with respect to the parameter 9. 

Finally the following assumption relates the pressure and the specific entropy defined below, 
through the variational principle. The variational principle is conjectured to hold for all models 
of statistical mechanics even if it is not proved in such a general setting. It is established for 
stable and finite-range interactions in Q and for infinite-range pair potentials that are regular and 
non-integrably divergent at the origin in [1, |^ . 

For any S G I and 9 G Q, the specific entropy of G is the limit 

hm (18) 

ra^-oo |A„| 

where Ia„ (P'^’®) = f In f^^{uj\^)P^^{duj) is the finite volume entropy of P'^’®. Note that this limit 
always exists, see [lOj. 

[VarPrinj: For any S G I and any 9 G Q, the pressure defined as the following limit exists and is 
finite 

p{5,9)= lim (19) 

n->oo |A„| 

Moreover the map 5 i—>■ p{S,9) is right continuous on I, i.e. lim5/_>5^5/>5p((5', 0) = p(5,9). In 
addition, for any 5, S' in I with S' > S, for any 9, 9' in 0, for all P^ G the inequality 

p{5,9) > -I(P'^'’®') - H\P^'’^') (20) 

holds. If {S', 9') = {S*,9*), then equality holds in (001) if and only if {S,9) = (S*,9*). 

The formulation of the variational principle presented here is slightly weaker than the classical 
one in statistical mechanics where the inequality (1201) is required for any probability measure P on 
CIt, and not only for P^ G ■ Note that the right continuity of the pressure holds trivially 
when S* is known, i.e. I = {<5*}. Let us also remark that the last part in [VarPrin] concerning 
(<5*, 9*) is just an identifiability assumption. 

We are now in position to state the consistency of the MLE. 

Theorem 1 Under the assumptions [Existence] , [Argmarx], [Stability], [MeanEnergy], [Bound¬ 
ary], [Regularity] and [VarPrin], for any {S*, 9*) G I xK. and any P G G^ , the MLE {5n,9n) 
is well-defined and converges P-almost surely to (S*,9*) when n goes to infinity. 

Example 1 (continued): The Strauss and the hardcore Strauss models satisfy all the assumptions 
of Theorem [T] which ensures the consistency of the MLE of S*, (3* and R*. The assumptions [Exis¬ 
tence] and [Argm8Lx] have already been discussed earlier and the stability assumption is obvious 
here since (j) > 0. The assumptions [MeanEnergy] and [VarPrin] involve a variational char¬ 
acterisation of infinite volume Gibbs measures which is standard in statistical mechanics [^. On 
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the contrary, [Regularity] is a technical assumption which is specific to the study of consistency 
of the MLE. Its proof is given in Section 15.31 when the hardcore parameter d* is positive and in 
Section when (5* = 0. Note that the Strauss model belongs to the large class of finite range 
pairwise interactions developed in details in Sections [4.11 and 


Example 2 (continued): The Lennard-Jones model also satisfies the assumptions of Theorem [T] 
The assumptions [Existence] and [Argmaix] have been cl^ked in Section 12.31 The stability 


assumption is already known, see for instance Section 3.2 in [25j, and is proved using the theory of 


positive type functions and Bochner’s theorem. The assumptions [MeanEnergy] and [VarPrin] 
are proved in a, a which deal with the variational principle in the general setting of infinite 
range pairwise interactions. Again the [Regularity] assumption is more specific and requires fine 
computations given in Section [S751 Note finally that the Lennard Jones model belongs to the class 
of infinite range pairwise interactions studied in Section 14.31 


As illustrated in the two examples above, all the assumptions of Theorem (TJ excepted [Reg¬ 
ularity], are natural and standard in statistical mechanics. They generally do not require much 
efforts to be checked since they are already proved for many models in the literature. In fact, only 
Assumption [Regularity] is really specific to the study of consistency of the MLE and needs a 
particular attention. Its proof strongly depends on the underlying model. Nonetheless, the tech¬ 
niques we develop in Sections [5.3115.41 and 15.51 in the cases of the Strauss, the hardcore Strauss and 
the Lennard Jones models (and their generalizations) provide tools and strategies to deal with it. 
We think that they could be easily extended to the study of other models. 

As a noticeable example stated in the following corollary, if the energy depends linearly 
on the parameter d, which is called the exponential model, then the assumptions [Argmax] and 
[Regularity] are automatically satisfied. Specifically, for any 1 < k < p, let (H^) be a family of 
interaction energies satisfying O- For any 0 S 0 and any bounded set A an exponential model 
takes the form 

Hi = ^9kHl ( 21 ) 

k=l 

Let us introduce the following modification of [MeanEnergy]. 

[MeanEnergy’]; The assumption [MeanEnergy] holds with Hq = where Hq comes 

from the decomposition ([T4|l applied to . 

Corollary 1 Assume that the family of interaction energies satisfy m and that 0 is an open 
subset of M2. Then under the assumptions [Existence], [Stability], [MeanEnergy’], [Bound¬ 
ary] and [VarPrin], for any {S*,9*) G I x )C and any P G , the MLE (Sn, 0n) is well-defined 
and converges P-almost surely to (6*, 9*) when n goes to infinity. 

Remark 1 In Theorem\^ there is no topological assumption on the set 0 except that 0^0. 
However, for identifiability reasons or in order to check [Regularity], some assumptions on 0 
will often come out, depending on the model. This is illustrated in Corollary [7] where we assume 
0 to be an open set. This assumption is usual and will still be used in the specific examples of the 
next section. Nevertheless this hypothesis is in general not necessary for consistency. Por some 
models, it is possible to check all assumptions of Theoreml^ when 0 is not open and 9* is on the 
boundary of Q. As an example, it is not difficult to prove that the estimation of (3* = 0 in the 
Strauss model (where J* = 0 and R* is known) is consistent, in which case 9 = j3, Q = [0,oo) and 
(I21|l is satisfied. This can be done as in the proof of Corollary QJ where the assumption Cll) in 
[Regularity] is verified by taking, for any rj > 0 and fio > 0, A/’(/3o) = {fio + 0/2} and r < 77 / 2 . 
Note that the same conclusion is not true if in this model we try to estimate both fi* and R* when 
(3* = 0, because then R* is not identifiable. Finally, the study of consistency for parameters on 
the boundary of 0 is specific to the model at hand and it is difficult to draw general results in this 
case. 
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4 Consistency for pairwise interactions and other examples 

4.1 Finite range piecewise continuous pair potentials with hardcore 

In this section, we illustrate how the general result of Theorem [T] applies to the class of finite-range 
pairwise interactions with hardcore. The following framework gathers all potentials of this type 
described in i, 0, 0 and [13 , including the hardcore Strauss model of Example I (detailed at 
the end of this section in Corollary [2]) , the hardcore piecewise constant pairwise interaction model 
and the Fiksel model. 

Specifically, we consider conditional densities given by 0 with a non-vanishing hardcore pa¬ 
rameter, i.e. i5min > 0, and a Hamiltonian defined as in ([2|) by 

hI{uj) = z Na{uj) + - y), (22) 

{x,v}euj,{x,y}nujA^I/) 

where z > 0, /I G and R = (i?i,... ,Rq) is a multidimensional interaction parameter in 
belonging to the subset TZ = {R £ R'?, 0 < i?i < • • • < Rq}. The pair potential interaction is 
defined piecewise as follows. For any x 

+ + if 1^1 Ri,R2,---,Rq, ,23) 

|min((/?f’^(|a;|), v?f)f^(|a:|)), if |x| = Rk for some k, 

0 R 

with the convention = 0. For any k = 1,..., o, we assume that there exists an open subset 
B of RP such that 

(/3, i?, x) !->■ ‘p^’^ilxD is continuous on B x TZ x R'^. (24) 

Note that from (1^51) and (|24p . the function {13, R) i—t (j)^’^{x) is lower semicontinuous which is 
crucial for the assumption [Argmax]. 

In this setting, the parameters to estimate are 5* and 9 = {z*,P*,R*) that we assume belong 
to / = [5min, oo) and 0 = R+ x B x TZ, respectively. For identifiability reasons, we further assume 
that 

Rl > 6* and X‘^ (x G R^ %.,+oo)(|a;|) (/*’^*(x) - /’^(x)) o) > 0 (25) 

for any {13, R) ^ (/3*,i?*). This implies that if {S,6) ^ (<5*,^*) then = 0 . 

Notice that the assumption Rl > S* is necessary to ensure identifiability of R* and S*, but we 
do not need to add the constraint i?i > i5 in the MLE optimisation ([5]) , as confirmed by the choice 
of K, in the theorem below. This is in particular important to agree with the setting of Lemma [l] 
where it is necessary that the optimisation ([8]) is carried out on a cartesian product / x /C. 

Theorem 2 Consider the Gibbs model defined by © with the Hamiltonian (EH) and the pair 
potential (l23l) . under the assumptions Jmin > 0, (l24l) and (l25l) . Let K be a compact subset of 
R+ X B xTZ such that {z*, {3*, R*) G 1C. Then the MLE of {5*, z*, (3*, R*) given by ([8]) is strongly 
consistent. 

Corollary 2 (Hardcore Strauss model) Consider the hardcore Strauss model of Example 1 
with parameters <5* > Jmin where (Jmin > 0, z* > 0, /3* G R and R* > S*. Let K be a compact 
subset o/R+ X R X R_|_ such that {z*, f3*, R*) G K. Then the MLE of {5*, z*, (3*, R*) given by ([8|) 
is strongly consistent. 

Proof. The hardcore Strauss model corresponds to the Hamiltonian ([2^ and the pair potential 
(1^ where q — p = I and = f3- In this setting TZ = R+. The assumptions and (051) thus 
hold trivially true for H = R and Theorem [5] applies. ■ 
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4.2 Finite range piecewise continuous pair potentials without hardcore 

We consider the same setting as in the previous section except that there is no hardcore, i.e. 5* = 0 
and S* is known, meaning that it has not to be fitted, or equivalently (imin = ^max = 0 implying 
/ = {0} in ([8|). The Hamiltonian it® is defined as in (1221) and (l23l) with the same continuity 
assumptions (1241) . The standard example we have in mind for this section is the Strauss model 
of Example 1 or any of the examples mentioned in the previous section without assuming the 
presence of a hardcore part. 

Whereas this setting may appear to be simpler than the previous one from a statistical point of 
view (there is one parameter less to fit), it requires additional assumptions to ensure the existence 
of the model and to control its regularity. This has been already pointed out for the Strauss 
model of Example 1 after the statement of [Existence] . Consequently, we further assume that the 
Hamiltonian iJ® satisfies [Stability] and that for any (/3, R) the pair potential is superstable, 
which means that it is the sum of a stable pair potential plus a non negative potential which is 
positive around zero (see (^1. The identifiability assumptions remain similar as in the previous 
section, that (1251) with d* = 0. 

To handle the second part of [Regularity] , we need a last hypothesis on the potential function. 
We assume that for any (/3o, Rq) G B x TZ and any rj > 0 there exists Pi G B{Po: ??) H H such that 
for any k = 1,... ,q 


^0uRo>^^^o,Ro_ ( 26 ) 

Theorem 3 Consider the Gibbs model defined by © with the Hamiltonian (1221) and the pair 
potential (1231) . under the assumptions ^min = ^max = 0, (1241) and (l2^ where S* = 0. Assume 
further that the pair potential is superstable and (1261) . Let K. be a compact subset o/ffi.+ xBxTZ 
such that {z*,P*^R*) G 1C. Then the MLE of {z*, P*^ R*) given by ([5]) is strongly consistent. 


This theorem can be easily 
hardcore) described in Q, [isj] . 


£plied to all standard finite range pairwise potentials (without 
20| and [27|. In particular, condition (l26l) turns out to be non- 
restrictive. Because of its central role in spatial statistics, we focus in the following corollary on 
the Strauss model of Example 1. 


Corollary 3 (Strauss model) Consider the Strauss model of Example 1 with parameters z* > 0, 
P* > 0, R* > 0 and there is no hardcore, i.e. 6* = 0 is known. Let K be a compact subset o/K^ 
such that {z*,P*,R*) G /C. Then the MLE of {z*, P*, R*) given by ([5]) is strongly consistent. 

Proof. Recall that the Strauss model corresponds in (l23l) to q = p = 1 and (pf’^ = p. As in 
the proof of Corollary [21 the assumptions ([2^ and ([25]) (where S* = 0) hold true with the choice 
B = M+. Moreover the associated pair potential is positive for any /3 > 0 and i? > 0, which 
shows that it is superstable and that [Stability] holds true. Finally the assumption ([26)) is verified 
if we choose Pi > Pq. Therefore Theorem [3] applies. ■ 


4.3 Infinite range pair potentials with a smooth parametrization 


In this section, we consider pairwise models in the spirit of the Lennard-Jones model given in 
Example 2. Our general setting concerns infinite range pair potentials without hardcore, i.e. 
<5min = i^max = 0, that are uniformly regular and non-integrably divergent at the origin in the sense 
of Ruelle [26| . Specifically, we assume that there exist two positive decreasing functions ip and y 
from R_|_ to K and tq > 0 with 


+ cx:) 


ip{t)C '^dt < 00, 


x(t)C '^dt = +00 
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such that for any parameter P £ B, where B is an open set in the function —>■ R 

satisfies for any a; G R'^ 

4>^{x) > xd^^l) whenever \x\ < tq and \(j)^{x)\ < t/’(|a;|) whenever |x| > rp. (27) 

The Hamiltonian iJ® is defined as in ([2|) by 

hI{u}) = z Na{u}) + ^ (t>^{x-y), (28) 

{a:,y}ei^.{a:,i/}nwA#0 

where 9 = {z, 13) is in 0 = R+ x B. 

We assume that (j)^ is a symmetric function on R‘^\{0} and that for any x £ R‘^\{0}, the map 
P !->■ <j)^{x) is differentiable on B. Denoting by \7(j)^{x) its gradient, we also assume that for any 
compact set 1C G B, for any P £ 1C and a: £ R"^ with |a:| > ro 

\V{x)\ <'il){\x\) (29) 

and , 

sup sup fI|x|<ro sup e~‘^^ max(^^(x), |V(/)^(a:)|) < oo. (30) 

p&Kxm'>- P'eic 

In addition we suppose that for any compact set 1C C B there exists an open set U in R^ and 
a stable pair potential p from R"^ to R such that for any u £ U, any x G R'^ and any P £ 1C 

Vp^{x).u > 4>{x), (31) 

where v.u denotes the scalar product of vectors u,u in R^’. 

Finally, for identifiability reasons, we assume that for any P ^ P' in B, 

^x G R'^, p^{x) ^ (x)^ > 0. (32) 

Theorem 4 Let {p^)p^B be a family of pair potentials which are uniformly regular, non-integrably 
divergent at the origin and satisfy assumptions (HID-®. Let 1C be a compact subset of 0 such 
that 9* = {z*,P*) belongs to 1C. Then the MLE of {z*,P*) given by ([8]) is strongly consistent. 

As a fundamental example, the following corollary focuses on the Lennard-Jones model of 
Example 2. Let the triangle domain T = {(n, m) G R^, d < m < n}. 

Corollary 4 (Lennard-Jones model) Consider the Lennard-Jones model of Example 2 given 
by (1281) and (jj) with parameters z > 0 and P = (A, B, n,m) £ B where B = R+ x R x T. Let K, 
be a compact subset o/R+ x B such that {z*,P*) £ K. Then the MLE of {z*,p*) given by (O is 
strongly consistent. 

Proof. Denote by ICj- the projection of K. onto T, i.e. (n, m) £ ICj- if and only if there exists 
z > 0, A > 0 and B such that {z,A,B,n,m) £ 1C. It is easy to find rg > 0 and two positive 
constants cg, ci such that pl^ given by (| 3 ]) is uniformly regular and non-integrably divergent at the 
origin for rg, 'p{t) = cot~"^° and x(t) = t~'^° with d < mo < m and d < ng < n for all (n, m) G /C 7 -. 
The map P 1 -^ p^ is clearly differentiable on B with 

V/(x) = (|x|-", |x|-™, -Aln(|x|)|x|-”, i?ln(|x|)|x|-™) . 

In adjusting the constant cg in the definition of p, we show that assumption holds. As¬ 
sumptions and are obvious. It remains to show ®. Let uq be the vector (1,1,1,1) in 
R"^ and u any vector in the open ball H(ug, 1/2) in R^. Denoting ICb the projection of 1C onto B, 
we find that for any P £ ICb and x G R"^ 

Vp^{x).u > i]I[o,i](|x|)(|xr" - |6gln(|x|)||xr"*) - ^lI(i.+oo)(|a;|)(ai + |6i|) IndxDIxr"* 
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where oi, bo and bi are chosen so that for any A and B, A < ai and bo < B < bi. There exists 
(5 > 0 such that for all (n, m) £ /C 7 -, m < n — S, whereby 

V(j)^{x).u > clI[o,i](|a;|)|a;|""'“ - c]I(i_+oo)(|a;|) ln(|x|)|a;|"™“ 


for some c > 0. The right hand term of this inequality is regular non-integrably divergent at the 
origin. Therefore it is stable and (1511) is proved. ■ 


We can also deduce from Theorem th e same kind of result as in 18| concerning exponential 
models. Note however that the proof in [l 8 | crucially relies on the convexity of /3 1 —>■ whereas 

the following corollary is obtained by different techniques. 


Corollary 5 Assume that in (1281) the family of pair potentials 'Is uniformly regular, non- 

integrably divergent at the origin, satisfies the identifiability condition dsn and that 


i=l 


where 4)i(x) > x(|a;|) if |a;| < ro while for i > 2, 4>i(x) = o{4)i(x)) when a; —>■ 0. Let K be a compact 
subset 0 /IR+ X B such that 9* = {z*,P*) belongs to 1C. Then the MLE of {z*,j3*) given by ([^ is 
strongly consistent. 

Proof. In view of Theorem SI we just need to check (l^ - dHTl) . Let /3 = (/3i,... ,/3p) G 1C and 
P' = {Pi, ..., Pp-i,Pp) £ 1C with Pp ^ Pp. Since is uniformly regular, we have \p^{x) — pl^ (x)| < 
2tp{\x\), meaning that \pp{x)\ < cV’(|x|) for some c > 0. We obtain likewise |</>i(x)| < ci/'dxl) for 
any i = 1,... ,p. Consequenlty \Vpl^{x)\ < ^ pc'I'(|x|) and (P^)) is proved. 

The relation (15(1)) holds trivially true by choosing rg sufficiently small. Finally, letting U be 
any open set included in 1C, we have Vp^{x).u = p'^{x) which is stable by assumption. H 


4.4 Examples of infinite-body interactions 

Our general result in Theorem [1] is also adapted to non-pairwise potential models. For instance, 
it is not difficult to generalize the results of the three previous sections to the case of a finite-body 
interaction of order greater than two, as for instance a triplet or quadruplet interaction. For brevity 
reasons, we do not include this generalization in the present paper. Instead, we give in this section 
examples arising from stochastic geometry which involve infinite-body interactions. Specifically, 
we focus on the area-interaction and the Quermass-interaction processes, but other models could 
have been considered as well. 

The area-interaction process [l[ , also called Widom-Rowlinson model in the statistical physics 
community (^ . is probably the most popular model of infinite-body interaction. For i? > 0, w £ If 
and A a bounded set in M.‘^, we introduce the notation 

.4«(cca) = A^ ( y B{x,R) 

\xGoja 

The Hamiltonian of the area-interaction process is defined for any a; € H and any bounded set 
A by 

Ha{^) = zN/^{uj) p (‘^A©B(0,2fi)) — -4^ (‘j-'A©B(0,2fl:)\A)) ) (33) 

where the parameter 6 = {z, R, P) belongs to 0 = M+ x ]R_|_ x M and the operator © stands for the 
Minkowski sum acting on the sets in Note that the unknown radius R is part of the parameters 
and is estimated consistently as stated in the following proposition. 

Proposition 1 LetTL = (FI®)a be the family of area-interaction energies defined in (l33)) and K, a 
compact subset of Q such that 0* = {z*, R*, P*) belongs to 1C. Then the MLE of 9* given by ([8|) is 
strongly consistent. 
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The Quermass-interaction model is a generalization of the area-interaction process, where not 
only the volume of the union of balls is involved in the Hamiltonian but also the other Minkowski 
functionals. We denote by fc = 1... d -|- 1, the d -I- 1 Minkowski functionals in and for short 

M^iujA)=Mk f U B{x,R) 

\x^UJA 

where R > 0, lu G and A is a bounded subset of R.'^. Recall that for d = 2, Mi corresponds to 
the Euler-Poincare characteristic, M 2 is the perimeter and M 3 is the area. We refer to Q for more 
details about Minkowski functionals. 

The Hamiltonian id® of the Quermass-interaction process is defined for any u; G H and any 
bounded set A by 



d+l 

Ha{uj) = zNa{uj) + ^ Pk (Aif (wA©B(0.2fl)) — Aif (wa©B(0,2B)\a)) (34) 

k=l 


where R > 0 , z > 0 and /3 = (/3i,..., Pd+i) G 

This model has been introduced in [15|. Its existence on R® has been solved so far only when 
d < 2 in 0. Therefore we restrict the following study to the case d < 2. Moreover, we assume 
that R is known and we only consider the MLE estimation of 9 := (z,/3) in (IMl) . The reason is 
that we did not succeed to prove that (flTl) in [Regularity] holds in presence of the Euler-Poincare 
characteristic Mi when R is part of the unknown parameter. With the assumption that R is 
known, the Quermass-interaction process becomes an infinite-body interactions exponential model 
and Corollary [T] applies. In this framework the set of parameters is 0 = R+ x R^. 

Proposition 2 Let R = {H^)a be the family of Quermass-interaction energies defined in (|34ll for 
d = 2 and let K, be a compact subset of Q such that 9* = (z*,/3*) belongs to K. Then the MLE of 
9* given by ([8]) is strongly consistent. 


5 Proofs 

5.1 Proof of Theorem [ 1 ] 

The proof is organised as follows. In a first step (Lemma [3]), we show that the hardcore parameter 
Sn converges to <5*. So, according to Lemma |2l it remains to prove that the minimizer of 0 —>■ 
) converges to 9* . This result is guaranteed via a general lemma on the convergence 
of minimizers of possibly non regular contrast functions (Lemma |4]). The final step of the proof 
therefore consists in checking the assumptions of this lemma, namely 

• the almost sure convergence of the contrast functions to a lower semicontinuous function 
admitting 9* as a minimizer (assumptions i), ii), iv) of Lemma |4]), 

• the control, when n goes to infinity, of the infimum of contrast functions evaluated on small 
balls with respect to the minimum of contrast functions evaluated only on a finite number 
of points (assumption u)). 

Note that Lemma [Hand the verification of its assumptions in the setting of Theorem |T] and for 
all models considered in Section |4] are the main contributions of the present paper. 

Let us now enter into the details of the proof. Let P be a probability measure in ’® and 
let w* be a realisation of P. If P is not ergodic, it can be represented as the mixture of ergodic 
stationary Gibbs measures, see [ 2 ^. Therefore the proof of the consistency of the MLE reduces to 
the case when P is ergodic, which is assumed henceforth. 

Let us start by proving the convergence of the hardcore parameter. 
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Lemma 3 The MLE 5n converges P-almost surely to 6* when n goes to infinity. 

Proof. From Lemma [U = min(5„, Jmax) where Sn is given by (flOl) . Clearly Sn is a decreasing 
sequence and (5„ > (5*. So it remains to prove that, for any S' > <5*, Sn is smaller than S' for n large 
enough. By the DLR equations ([S]) and the definition of local densities d?]) , there exists a bounded 
set A such that 

P(there exist x,y £ uja such that (5* < |a; — ?/| < 5') > 0. 

Thanks to the ergodic theorem, this implies that 

P(there exist n> 1 and x,y £ such that d* < |a; — ?/| < 5') = 1 
which proves the expected result. H 

According to Lemmas [5] and [31 we have now to prove the consistency of where 

On = {ujIJ. 

The strong consistency relies on a minimum contrast function result stated in the next lemma, 
where neither the contrast function nor its limit need to be continuous with respect to 9. 

Lemma 4 For any 6 £ Q, let (/i® )n>i be a family of parametric measurable functions from to 
R. Let K, be a compact subset of Q and let P be a probability measure on LIt- We assume that 

i) For any 9 £ Q, hn converges P-almost surely to a finite real number (denoted by h^) when 
n goes to infinity. 

ii) The function 9 h^ admits a unique minimum over 0; 

9* = argming^Qh^, 


and 9* £ 1C. 

Hi) For P-almost every u> £ fir and for n sufficiently large, 9 ^ hn admits an infimum over 
1C attained by at least one element 0„(a;). 

iv) The function 9 ^ h^ is lower semicontinuous on 0 (i.e. for any 9 £ Q, Wm.ra.igi^g h^ > 

h^). 

v) There exists a function go from M+ to R+ satisfying 

lim go{x) = 0 

and such that for any e > 0 and for any 9 £ 1C, there exists a finite subset M{9) included in 
B{9,go{e)) fl 0 and r{9) > 0 such that 


P 


lim sup 

n—¥oo 


( min 

9'GAr{9) 


h 


9' 

n 


inf 

9'£B{9,r{e))nK 




(35) 


Then the sequence = argming^f^hn{9) converges P-almost surely to 9* when n goes to 
infinity. 


Proof. Let e > 0 and be the following set 

Ae = £ /C © B{0, (7o(e)) such that h^ < h^ + 2e| , 

where go comes from v) and B{0,go{s)) denotes the closure of the ball B{0,go{s)). By iv) the 
function 0 i—>■ /i® is lower semicontinuous and so the set is a compact set. Moreover by ii), 


{ 0 *}= 

£>0 


(36) 
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We deduce that the diameter of goes to zero when e goes to zero. Let us consider the set 
Be := Ag © B{0,go{e) + e). Clearly its diameter tends to zero as well. So to prove that (0„) 
converges P-almost surely to 9*, it is sufficient to show that for any e > 0 

P (limsup e ) = 0, (37) 

where PJ = /C \ (Be H /C). 

Since K. is compact, there exists a finite sequence (0i)i<i<N in K. and {ri}i<i<N in 1R+, where 
Ti = r{9i) comes from v), such that the union of balls B{9i,ri) covers JC. Note that if (1351) holds 
for some r, then it holds for any r' < r. So without loss of generality we can assume ri < e for all 
i = 1,... ,N. Let Ie be the subset of indexes i G N} such that the ball B{9i, ri) intersects 

the set Bg. We denote by A/) the subset Af{9i) in assumption v) and by the union of A/) for 
i G le- By definition and since for alH = 1,..., A^, < e and A/) C B{9i, go{e)), we have 


B^ C [J B{9i,ri), A4 C/C © B(0,go(e)) and A4 n Ae = 

i 6 le 

The two latest properties imply that 

V6»gA4, 

Let us finally consider the following sequence of random variables: 


6^ = max ( min ht — 


ieXe \6»'e7Vi 

Note that from (155)) and since is finite 


inf / 
S'GB{9i,ri)nK 




(lim sup > ej ) = 0 . 

\ n^oo ^ ^ J 


To prove ()37l) . we consider the following inclusions, 
lim sup < G Bg i C lim sup < inf 

n—¥QO ^ ^ n—¥QO [ 


C lim sup < min /i® — /i® < <5^ 
n—¥cc) 1 eeA4 


C 


li^sup {(5^ > U -hn <e 

= lim sup < <5^ > e M I lim sup < min < e 

n—^oo ^ ^ n—¥oo e&Ne 


c 


lim sup < (55 > e M > lim inf ( min — h', 
n^oo I J ^ I yeeATe 

lim sup |(5® > M i min < e| M ^ 

n—^oo t J ^ (6»e7Ve J ^ 


(38) 


(39) 


(40) 


where J^e is the event where /i® does not converge to , when n goes to infinity, for some 9 in 
Afe U {0*}. Since A4 U {0*} is countable and by assumption i), P{^e) = 0. Moreover, from (1551) . 
the set {mineejN/j < e} is empty. Recalling (HOI) , we thereby deduce (1571) and the lemma 

is proved. ■ 

Let us show that the assumptions i), ii), in), iv) and v) of Lemma |4] hold for the family of 
contrast functions 

Hztf) . KMkJ 




|A„ 


|A„ 
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To prove i), first note, from [MeanEnergy], [Boundary] and the ergodic theorem, that for 
any 0 £ 0 and P-almost all ui 

hm = Ep{H^,) = H\P). 

n-^oo |A„| 

Second, recall from assumption [VarPrin] that for any S G I and 0 £ 0, ln(Z^’^)/|A„| converges 
to p(S, 0) and S i—^ p(S, 0) is right continuous. On the other hand, from Lemma |3l 6n —^ almost 
surely with i5„ > S*. Let 0 £ 0, e > 0 and 5+ > S* such that \p{S*,0) — p(5+,0)| < e. Let n be 
sufficiently large so that (5„ < <5+, | ln(Z^^’®)/|A„|—p((5+,0)| < e and | ln(Z^J®)/|A„|—p(5*,0)| < e. 
Then, since 6 i—>■ is decreasing, we have 


ln(zi-^)/|A„| -p(r,0) < ln(zf/)/|A„| -p(5*,0) < 


and 


p{S*,e) - ln(Z^";")/|A„| < pi6*,e)-p{6+,0) +p(5+,0) - ln(zf’®)/|A„| < 2e 


proving that ln(Z^”’^)/|A„| converges to p{6*,9) almost surely. 
Hence for any 0 £ 0 and for P-almost every w* 


lim = h 


where 

=p{S*, 9)+H^{P) 

which proves assumption i). 

To prove Assumption ii), note from [VarPrin] that for any 0 £ 0, any P^ and any 

p5* ,9* ^ qS* ,9* 

p{5* , 0) > -I(P^* ) - P® {P^" ) 

while -I(P^*’®*) = p{S*,9*) + so that 

/i® > h^\ 


From [VarPrin], the equality holds if and only if 0 = 0* which proves Assumption ii). 
Assumption Hi) is given by Lemma |2J 

The function 0 i— H^{P) is continuous thanks to (fTOl) in assumption [Regularity]. By as¬ 
sumption [VarPrin] we have for any 0 and 0' in 0 


p{S*,9') > p(<5*,0) -f P®(P‘^*’®) -P®'(P'^*’®), 


where P^ G ■ By continuity of 0 i—?► H^{P), it follows that 0 i—7> p{S*, 9) is lower semicontin- 
uous. Therefore 0 i—^ /i® is lower semicontinuous and assumption iv) holds. 

It remains to prove assumption v). Let us start with some preliminary results. 

Thanks to [MeanEnergy], [Boundary], (TlBl) in [Regularity] and the ergodic Theorem, we 
have that, for any 0 £ 0, any r > 0 and P-almost every w. 


lim sup T-—r sup 
rn-oo |A„| 0i^f^ 

\e-e'\<r 


HtiuJAj 




< 9{r)- 


(41) 


This inequality allows us to control the variation of the infimum of the specific energy in the 
contrast functions. The following lemma deals with the variation of the finite volume pressure. 
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Lemma 5 For any rj > 0 and any 6 q G K, there exists a finite subset Af{0o) <Z B{0o,r]) D Q and 
0 < r{0o) < rj such that for n> 1 and any S € I 


1 


ln(Z^’^) — inf 


1 


ln(^Af) < e'"+i5(r(6»o)), 


(42) 


eaN^eo) |A„| eG-B(eo,r(eo))nK: |A„ 

where g comes from assumption [Regularity] and k from assumption [Stability]. 

Proof. For any S G I and any 0, O' in 0 

ln(Z^’® ) - In(Z^f) = In 

= In (^y e-< (-A„) (C.A(dwA„)) 

= In Eps^e (e-^A„ +Hi ^) ( 43 ) 

where is the probability measure (1/Z^’^)e~^''n Hqs t^a„- By Jensen’s inequality we have that 

ln(4f') - ln(4f) > Ep^AHl - H^J. 

For any p> Q and 0o in Ai we choose Af{0o) and r{0o) > 0 as in ([T7|) . Denoting [x]+ = max(x, 0) 
and writing for short M for J\f{0o) and B for B{0o,r{0o)) r\ JC, we obtain, for any 6 > Jmin, 

min ——- \n{Z^ff) — inf ——- ln(Z^’®) 
eeA^|A„| ^ eeB|A„| ^ 


= mm . 

Q'^B |A' 




- - Hij) 

deAfg'^B |Aral A A„ / 


|A 


1 


, , , min sup E^s.e 
|A„| eaM g,^B 




■Na„ 


< T -—- min sup sup 
|A„| e'GBuiAn^^^L 


NAr. 

hUujaJ - HUujaJ 


^ (maxSpi.s (A^a„) ) 

|A„| \eeAr "7 


1 


( 


maxEpS,e {Na„) 


|A„| \e^Af ^a! 


A^a„ (wa„ 


mm sup sup 
eeAf e'eBcc;A„ea^ 


— max inf inf 
deAf ®'6Bi^A„6n^ 


EpS,e {Na„) 


PaA^AJ - gA„KJ 

A^a„(wa„) 

( HUcoaJ - H^^JcoaS 

\ Na„ (wa„ ) 


< 


y p„, 

|A„| \B&Ar 


maxEp6,e {Na„) g{r{0o)). 


(44) 


Let us control EpS,e{NA^) by entropy inequalities. By definition of the entropy, the assumption 

■A-n 

[Stability] and the standard inequality > 7rA„({0}) = we find that 


^aAp'a!) = I In 


J A 


7S,e 


<-ln(4f)-y <„Pt' 

^ |An| + Kpp«,6i (A^A„)- 


(45) 
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The entropy is also characterised by 


(/ 9Phi 

where the supremum is over all bounded measurable functions g, see for example [ 2 ^. Choosing 
g = {k + 1)Na^^Na„<C for any constant C > 0 , we get 

IaApI:!) > (« + 1)Ep^ANa„1na„<c) - In (^J 

> {K + l)Ep^,e{NA„lLNA^<c) - (e"+i - i)|A„|. (46) 

Let C goes to infinity and combine (HSl) and (HSl) . we obtain that for any {S,9) 

EpS,e{NAj<e^+^\Ar,\. (47) 

An 

The inequalities (l44ll and (lT7)) imply (l42ll . ■ 

We are now in position to prove v). The function go in v) depends on g in dU) and Lemma [S] 
as follows. For any x > 0, we choose go{x) > 0 such that 

5 ( 250 ( 3 ;)) + e'^+^5(go(a;)) < x. 

This choice is always possible since g{u) tends to 0 when u goes to 0. 

Let £ > 0 and 9q G K. and consider the finite subset Af{9o) C B{9Q,r]) fl 0 and the positive 
number r{9o), with r{9o) < 77 , given by Lemma [S] where 77 = go{e). In the following, we write for 
short Af for A/’(5o) and B for B{9o,r{9o)) fl/C. Using (ITO . Lemma[5]and the definition of go, we 
have that for P-almost every w, 



limsup ( min /7®(a;A„) - inf h®(a;A„) ) 

n—yoo J 

“ “AfS”ra (""-<"A.) + UZi-f)) - M {him + ln(Zt-“))) 

- <“A.) - “ "L(“a.) + mnlKZf:'") - “ 


< max limsup 

dGhf n—>00 


h'-nl 9'eB 


HijuJAj 


HiS^hJ 


■ lim sup 


1 

1^ 


( minln(zf"’^) 


inf ln(Zf"’^) 


< 5 ( 2 5o(£)) + e'"+^g(5o(e)) < e, 


which proves v). 


5.2 Proof of Corollary [T] 

To apply Theorem [U we have to prove that [Argmarx] and [Regularity] hold true in the setting of 
exponential models. Assumptions [Argmax] and the inequality (fTHll in [Regularity] are obviously 
satisfied since the energy 77^ and the mean energy iJg are linear in 9. 

It remains to check dm in [Regularity]. Let /C be a compact subset of 0 and e > 0 such 
that /C 0 i?(0, e) C 0. For any 77 > 0 and 9o G 1C, we choose r sufficiently small such that 
5 := (1 + r/e)9o G B{9o,g) D 0 and we set Af{9o) — {5}. Note that for any 9' G B{9o,r), 
f:{9 — 9') G B{9o,t) C /C © B{0,e). Then from (PT|) and the assumption [Stability], for any 
9' G B{9o,r), for any wa G 

^{HiicuA) - h({uja)) = nf-^'P^ujA) > -i^Na{u,a) 
which implies (77® (wa) — 77® {(jJa)) > — rf 7Va(wa) and thus (fT71) . 
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5.3 Proof of Theorem [2] 

We check the assumptions of the more general Theorem [T] with the choice Q — Z x T where Z and 
T are bounded open subsets of (0,oo) and B x TZ respectively satisfying Z x T D 1C. We denote 
in the following 6 = {z, /3, R), where z C Z and (/3, R) G T. 

For any 6 £ I and (/3, R) £ T, the pairwise interaction ooll[o^ 5 ) + is superstable, finite 
range and bounded from above. The existence of Gibbs measures in this setting, i.e. assumption 
[Existence], is proved for example in [^. Moreover, we can deduce the following useful result 
from Corollary 2.9 in [2^: for any P £ , for any bounded set A and any c > 0, 

Ep < 00. (48) 

Superstability implies [Stability] , where the lower bound can be chosen uniformly over 0 since 
this set is bounded and using the hardcore property. The decomposition in [MeanEnergy] can 
be done by 

hL (^) = E ^0 o (^) (49) 

where 

Hq{uj) = zN^^(w)+ ^ + ^ (50) 

{a:,y}eu;Ao 

and 

The finiteness of the mean energy given by (US is easily deduced from the finite range and 
boundedness of along with (l4^ that implies finite second order moments for P. On the other 
hand from (I2|) we have 


x€UJAn 

Therefore [MeanEnergy] and [Boundary] reduce to prove that for P-almost every w £ Clp 

1 


sup 
iP,R)&r 


E 


\x - y) 0. 


(51) 


xGuJAn ^y^^A^ 


Recall the definition of In = {—n,... ,n — 1}'^. Since has finite range and is continuous, 
there exist two constant c and f large enough such that 


1 


IA _ 

(/ 3 ,fl)eT I n| 




ieln\ln-l 


< 


\K 




ieln\ln-l 


From (1481) . for any P £ Q^’^, Ep{Ng^. -j) is finite and by stationarity Ep{Ng^. -^) = Ep{Ng^Q p^). 
The ergodic theorem thus implies that P almost surely 


iXl E ^Pi^B{0,f)) 

I ”1 *e/„ 


which yields (ED). 
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Assumption [Argmax] holds because (/3, R) i—>■ is lower semicontinuous by construction 

and the number of discontinuities are finite, so it is easy to change the sense of the brackets in (051) 
to obtain an upper semicontinuous version. 

It remains to check (a) [Regularity] and (b) [VarPrin]. 

(a) We start to prove (flBl) in [Regularity]. Note from (l50]) that 

Hq (w) - (w) = {z- z) Aao (w) 

{x,y}ei^AQ a;ea;Ao 

If \9 — 6'\ < r, obviously there exists a positive function g with lims_,.o g{s) = 0 such that 

\z- z'\Nao{uj) < g{r)NAo{uj). (53) 

The two sums in the right hand side of (1551) can be handled similarly. We give the details for the 
first one only, which reduces from (1231) to consider the following generic term for P-almost every uj 

X \'^iRk-i,Rk}(\^ - y\) - y\) - ^{R'^_„R'^)i\^-y\)y^k i\x-y\)\ 

< X ’■” l(l^-yl)+c X 

{x.ylgWAo {x,y}Goj^g 


where c is a positive constant that comes from the boundedness of Since the map (/3, R, x) !—>■ 
^k’^{\x\) is continuous it is uniformly continuous on any compact set and there exists a positive 
function g, independent of x, y and 0, such that lims_>o g{s) = 0, and if |0 — 0'| < r then 

X \^k^-^k'^\‘^x-y\) ^ X 5(’’) < 5(O^Ao(w)- (54) 

{a:,i/}gWAQ {a:,i/}et^Ao 


On the other hand, denoting Ai?^ the symmetric difference (i?fc_i, i?fc)A(i?^_^, iij.), we have 
for any R € TZ 



( \ 


( \ 

Ep 

sup X ]lAflfc(|x-y|) 

= Ep 

X X lAijJlx-yj) 






( 


— Ep 


c^-h {x\uj) 


' Ao 


X 


sup 


llAflfcdx - y\) dx 


where the last equality comes from the Georgii-Nguyen-Zessin (GNZ) formula UM and h^{x\uj) 
denotes the local energy needed to insert the point x into the configuration w, defined for any 
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A 3 a; by h^{x\uj) = H^{u} U a;) — Applying again the GNZ formula we obtain 

/ 

Ep sup E lAflfc(|a:-2/|) 


/ 


= Ep 


[ sup - y\)dxdy 

JA?. R'gTZ 


R'&TZ 

\R-R'\<r 


< i?P (e'^^^offisco.s')) f sup yyARk{\x - y\)dxdy 

JAg K'eTC 

< Ep (e‘=^'^oeB(o,*')) IAoIIAWI 


(55) 


where c and f are two constants large enough and Ak(r) is the union of the two rings {y € 
R'^, Rk-i — r < \y\ < Rk-i + r} and {y, Rk — r < \y\ < Rk + r}. Thanks to (HSl) and since 
|-4/c(r)| —>• 0 as r —^ 0, the upper bound in (1551) tends to 0 as r —>• 0. The same conclusion holds 
true in (1551) and (1551) because Ep{Nao) < oo and Ep{N\^) < oo. The combination of these results 
shows (ITbl) . 

Let us now prove (ITTl) in [Regularity], Let rj > 0 and 9o = {zo,Po,Ro) € /C where Rq = 
{Ro.i, ■ ■ ■, Ro,q)- We choose r sufficiently small to ensure that 

(i) r <T], 


(ii) B{eo,r) C 0, 


(hi) 2r < inffe=i,,,,,q |i?o,fe - -Ro,fc-i|, 

(iv) for all fc = 1,..., g, if (Ro,fc) 7^ i^o,k), then for all 9', 9" G B{9q, r) and for all x 

such that ||x| — i?o,/c| < f, we have 

(|a^l) - (|a:|)) {fk°’^°i^o,k) - ^t^i°{Ro,k)) > 0, 

which is possible by continuity of {^,R,x) !->■ ip^’^{\x\). 

Then we fix Af{9o) = {{zo, Po, R)} where R = (i?i,..., Rq) is defined as follows : if (/7^°’^“(i?o,fc) > 
ip^°^^°iRo,k) then Rk = i?o,fe + r, if ip^°'^°{RoP < (/?f+’f“(i?o,fe) then Rk = i?o,fc - r, and if 

^f°iRo.k) = ^t'^+i%Ro.k) then Rk = Ro,k. 

For 9 = {zo,Po,R) and 9' G B{9o,r), the proof of (llTl) amounts to control 


HjiiXA) - <(^ a ) 
Aa(wa) 


{zq — z') + 


1 

Na{u}a) 


E 

{x,y}^uj,{x,y}nuiA^^ 



y) - '^\x - y)'^ . (56) 


The first term (zq ~ z') is greater than —g{r) for a positive function g with lims^o5(s) = 0. For 
the second term, let us introduce the middle points Ro.fc = ^(^o,fc-i + Ro,k) for fc = 1,..., g — 1, 
,Ro,o = 0 and Rg.g = Ro,q + 2r. Note that from our choice of r, for any k both Rk and belong 
to [,Ro,fej-Ro,fe+i]- From (l23l) . we can write 


^Po.R^x) - • 


k^O 

= E {^lRo,k,R<c)i\^\'l‘y’k°"^i\x\) + I[fi..flo,.+i](|a;|)<pf+f(kl) 

k^O 


( 57 ) 
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If (Ro^k) > then Rk > R'f. and the A:-th term in the sum above writes 


- V>1 (|a:|)) + (v5f“’^(k|) - <P^k+i (Ia^l)) 

+ (|a;|) - (|a:|)) . 


From our choice of r, see (iv), the second term is always positive, and by uniform continuity of 
x) I—>■ (p^’^{\x\), there exists a positive function g, independent of x, with lims_,.o g(s) = 0, 
such that the two remaining terms are greater than —g(r)lIjQ _^](|x|). We thus obtain that the 

fc-th term in ((Ffl) is greater than —g](l^l)- 


If {Ro,k) < ^^+ 1 ° (^o.fc)) we obtain the same lower bound by using the fact that Rk < R'^,- 

If {Ro^k) = then by uniform continuity ^k^{\x\) - (|?/|) < g{r) for any 


k' = k^k + 1, any |x|, |j/| in [-Ro.fej -Ro,fe+i] and any d, 6' in B{0Q,r), and the same lower bound holds 
for the the fc-th term in (157)1 . 

Coming back to (1551) , we deduce that 


HUua) - <( cca ) 
7Va(wa) 


> -gir) - (g+ 1 ) 5(0 


1 

Na{uja) 


{x,y}eu:,{x,y}nuJA^<ll 


y\)- 


For any oj G there exists c > 0 such that for any x G w, yj^x lI[o, Ro,q] (|x - 5 I) < c 

because of the hardcore distance dmin > 0. Consequently the lower bound above is greater than 
—g{r), up to a positive constant, and this completes the proof. 

(b) The variational principle in this setting is proved in , which implies (1201) with equality if 
and only if P^'. When (d',5') = {5*, 9*), P‘5*’®* = P^’® is equivalent to {5,9) = {5*, 9*) 
from our identifiability assumption (I25F Hence, in order to verify [VarPrin], it remains to prove 
that the pressure ca is right continuous in 5. 

For any d > 0 and 5 G 0, from the rescaling oj ^ Soj, 


r^s,e 

^An 


= g-|An|(l-l/5 ) J 





where = / e {uj a)'^ k{doj a) is the partition function associated to the Gibbs mea¬ 

sure defined on with energy function H^^{uja) = H^AiH^A)) ~ 1Va(w) ln(d). Denoting p{S, 9) 
the pressure associated to the latter measure, we deduce 

p{6,9) = ^p{5,9)-l + l/d‘^. (58) 

The energy function follows a decomposition like ((TT)) where 

Ho\uj) = {z- ln((5)) Nao{uj) + ^ /’®*((5x - Sy) + ^ ^ cj)^’^{5x - Sy). 

{x,y}GoJSAo xe(XSAQ,y£‘^5l\^ 

By similar arguments as used earlier in (a) to prove (ITbl) . we deduce that d i-!> Ep{Hq^) is contin¬ 
uous. Moreover, inherits the superstable, lower regular and regular properties of 0^’^, which 
shows from that a variational principle holds, i.e. for any {5,9) and {5',9') 

p{5',9')>p{5,9) + H^'%P^'^)-H^''^'{P^'^), 
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where denotes a Gibbs measure associated to . Note that contrary to (EUl) . this inequality 
is valid for any S,6' because P^’^ is supported on for any {6,0). We therefore deduce that 
(5 I—>■ p{S, 9) is lower semicontinuous, and so is <5 1 —>■ p{6, 9) from (l58l) . Since 6 i—>■ is decreasing, 
6 1 —^ p{6,9) turns out to be both decreasing and lower semicontinuous, which shows that it is right 
continuous. 


5.4 Proof of Theorem [3] 

As in the previous section, it is sufficient to check the assumption of the main Theorem [TJ The 
assumption [Stability] is included in the assumptions of Theorem|3]and has not to be showed. The 
proofs of assumptions [Existence], [MeanEnergy], [Boundary], [Argm 2 Lx] and the first part 
of [Regularity], i.e. (fT6l) . are the same as in the proof of Theorem!^ The variational principle is 
shown in for any finite range stable pair potential. Since / = {0}, the right continuity of the 
pressure is not required here, so Assumption [VarPrin] follows. It remains to prove the second 
part of [Regularity], i.e. (1171) . 

Let r? > 0 and 0o = (zojPo, Ro) S fC where Rq = {Ro,i, ■ ■ ■, Ro,q)- We denote by R the maximal 
value of Rq when 9 = {z, /3, R) G JC. We choose r > 0 small enough such that 


(i) r < 77 , 

(ii) B{0o,r) C 0, 

(hi) 2r < inffe=i,...,, \Ro,k - Ro,k-i\, 

(iv) there exists Pi € B such that for any 9' G B{9o,r), any 9 = {z,P,R) G B{{zo, Pi, Rq), y/qr), 
any k = 1,... ,q and any x G B{0, R) 

pi^'^\\x\) > ’■”^(1x1), 

which is possible by continuity of {P,R,x) !->■ p^'^{\x\) and thanks to the local positivity 
assumption (l26l) . 

(v) forallfc = 1,... ,q,iiipl°’^°{Ro^k) p (-Ro.fe), then for all 6»', 6(" G B{9Q,r)VJB{{z(i, Pi, Rq), 

and for all x such that ||x| — i?o,fe| < r, we have 

(|a;|) - V>lPi (|a:|)) {Ro,k) - <^fe+f“(Ro.fc)) > 0, 

which is possible by continuity of {P,R,x) !->■ ip^'^{\x\) and since Pi can be chosen as close 
as we want to Pq. 

(vi) for all k = l,...,q, if = ^l+i°{Ro,k), then for all 6»' G B{9o,r), 9 G B{{zo, Pi, Rq), y/^) 

and for all x such that ||a;| — i?o.fe| < r. 


{P,R) 


{\x\)>‘pi+f\\x\) and (kl) > ’"^(|a;|). 


d/3'.fl') 


which is possible for the same reasons as in (iv). 


Now we choose A/’(0o) = {{zq, Pi, R)} where R = (i?i, i? 2 , ■ ■ ■ ,Rq) is as in the proof of Theorem 
[21 see (l56l) and before. Note that {zq,Pi,R) is in B{{zo, Pi, Rq), y/qr) and that Pi can be chosen 
sufficiently close to Po to ensure J\f{9o) C B{9Q,r]). Following the same calculus as in ([57|) . we 
obtain that for any 9' = {z',P',R') G B{0o,r) 

0 /3 iA(2 ,) > (j)^ {x). 

We deduce that 
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and the second part of [Regularity] is proved. 

5.5 Proof of Theorem [5] 

We check the assumptions of Theorem [TJ For any P G B the pair potential (j)^ is regular and non- 
integrably divergent at the origin. As it is well known (see [251 ] for instance), (j)^ is superstable and 
the existence of an associated Gibbs measures follows, i.e. [Existence] holds true. Superstability 
implies stability which shows assumption [Stability] . Since the map {3 ^ is differentiable and 
so continuous, assumption [Argmaix] is obvious. The decomposition in [MeanEnergy] is done 
as in In particular the mean energy (ITSl) is 


H%P) 


Ep 


Nao + 


{x,y}eui^g 



(59) 


where 9 = {z,j3). To check assumptions [MeanEnergy] and [Boundary], we thus have to prove 
that (l5^ is finite and that for any compact set K, d B, 


sup 


I A 




0 . 


(60) 


For this purpose, we need the following Ruelle estimates. 

Proposition 3 (|26l]I We define, for any i G ijji = 'if(max{\i\ — 1,0), where the function 
comes from assumption (]27|) . Then < +oo and for any c > 0 


Ep ^e‘^'E,iezd ^ + 00 . 


Let us show (l60ll. Note that 


sup 


^ I A 


E {x - y) < + A 2 


X^UJAn 


with 


E sup|/|(a;-y) and A 2 =-E ^ sup |/|(a; - y), 

y^uj,y^x 

\x-v\<ro \x-v\>'ro 


where ro comes from assumption (EZj) and no is an integer greater than rp. 
By the spatial ergodic theorem, P-almost surely 


lim 

n—^oo 


1 

1^ 


V sup - y) = E 

ye“. 

\^-v\<’ro 


E sup|<()^|( 

\x-y\<ro 
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Therefore, Ai goes P-almost surely to 0 if the expectation above is finite. By the GNZ equation, 
stationarity of P and assumption (EZl) 


E 


V snp\(l)>^\{x-y) 

\x-y\<ro 


V sup/(2/) 
V ye^ /3e/c 

|y|<’'o 


< z*e( e sup (/)^(?/)'] 

V y^u, /^e'C J 

|y|<'ro 

< sup ) 

V ^ PeK J 

\v\<i~o 

< z*e( sup sup sup 

V ’ J pe>Cy^B{ 0 ,ro) P'eK 


(61) 


which is finite by assumption (1301) and Proposition [3] 
Regarding the term A 2 , for any integer K we have 


A2 ^ \ , I '*/’|z-j|-^Ti(Ao)(^)-^Tj(Ao)(^) 

I A I~ \ A r- \C 




< 


|A„ 


A^Ti(Ao) (^)A^Tj(Ao) (^) + 1 ^ I V'li-jl A^Ti(Ao)(w)A^rj(Ao)(^)- 


ieA„\A„-K,. 


leAn.K.jez" 

|i-*l>x 


(62) 


By Proposition |3l A^Aq X)jez<i '4’j^Tj(Ao) i® integrable which implies thanks to the ergodic theorem 
that, when n goes to infinity, the first term in (1531) goes to zero and the second term goes to 


E 


Nao 2^ V’iA^T,(Ao) 
jeZ'i,Ul>K 


Choosing K large enough this last expectation can be made smaller than any positive value. 
Therefore the term A 2 tends P-almost surely to zero as well and (1501) holds true. 

By similar computations as above, using the Ruelle estimates in Proposition |3l we obtain that 
the mean energy (l59l) is finite. Assumptions [MeanEnergy] and [Boundary] are verified. 

The existence and finiteness of the pressure in assumption [VarPrin] is proved in Theorem 0.2 
of [ 2 ^ . The variational principle for such Gibbs measures is proved in [^, and thanks to the 
identifiability assumption (1301) we get [VarPrin]. 

It remains to prove [Regularity]. Let us show (TTOl) . For any compact set Z x 1C C & with 
Z C R+ and K. C B, for any r > 0, 

< rEp{NAo) + Ep { E E sup — \{x — y) 

< tEp^Nao) +rEp ( E E sup \V(j)^\{x — y) 

a;ecJAo yeci;\{3;} 

The proof of (1161) is completed if we show that the two expectations above are finite. This is clear 
for the first one. We split the second one according to |?/| < ro or jyl > tq where xq comes from 


Ep 


sup 

e'&zxK 

\9-B'\<r 


TT& TTU 

Iln — -tin 
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assumption (EZl). Thanks to similar computations as in (EB) we get 


Epi sup\V(j)^\{x-y)] 

|y|<’'o 

< ro) ] sup sup sup | | (j/) 

V ’ J /3elCy&BiO,ro)y'eK 

which is finite by (l30l) and Proposition [U while from (l2^ 

M E E sup \V(j)^\{x - y)\ < e(n^„ E V'jiVr,(Ao)) 

^ a:ec-;Ao yec^\{x} ^ ^ ^ 

\v\>xo 

which is also finite. This proves (ITSl) in [Regularity]. 

To prove (113, let ?7 > 0 and Qq = (zo,/3o) G 2 x /C. There exists 0 < ^ < y such that 
E. Consider the open set U associated to /C©R(0,5) through the assumption (IHTl) . Let 
us fix uo G and e > 0 such that B{uo, e) is a ball included in U. We choose r sufficiently small 
to ensure that r < ^ and ^ := /3o + f mq belongs to B{l3o, y)- By assumption, for any /3' G R(/3o, f) 
and x G there exists t G [0,1] such that 

(j)^{x) — (j)^ {x) = +‘(/3-/3 )(x).(/3 — j3'). 

Note that the above choices ensure that jS' + t{j3 — j3') belongs to i?(/3o, C /C © R(0, and that 
^(/3 — j3') belongs to B{uo, e) C U. Hence by assumption (|5T]) . for any x G 

-{(j)^(x) - (j)E{x)) > ^{x) 
r 

where ^ is a stable pair potential. We finally choose A/’(0o) = {0} := {(zo,/3)}- Therefore for any 
bounded set A, any 9' = {z', 13') G [zq — r, zq + r] x B{l3o,r) and any configuration w, denoting by 
A > 0 the stability constant of we have 

hI{u}a) - (wa) > (zo - z')Na{u}a) + ^ - y) >-r {I + A/e)NA{ujA) 

{x,v}^u>x 


which proves (El) and concludes the proof. 

5.6 Proof of Proposition [1] 

We apply Theorem [T] The local energy h^{x\oj) associated to the area-interaction process is, for 
any A containing x, 

h%x\uj) = Hiiuj Ux)- Hliuj) = z + /3AN B{x, R) \ |J B{y, R) j . (63) 

Since /C is compact, z = sup^g^s^ z, R = snpg^j^R and j3 = supgg^g |/3| are all finite and hP is 
uniformly bounded: 

sup sup sup \h^{x\uj)\ < z + ^(2R)'^. 
xgR'^ ee/CcjGn 

This local stability property implies [Stability]. Since h® has also finite range 2R, an associated 
Gibbs measure exists and the variational principle holds with the existence of the pressure, see 
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24 1 and [^. This gives [Existence] and [VarPrinj. Since the application 6 i—is con¬ 
tinuous, the assumption [Argmax] is obviously satisfied. As for [MeanEnergy] we suggest this 
decomposition (other one could have been considered) 




with 


= zNa,{uj) + f Ao n U B{x, R) j . 

\ X^CJ / 


(64) 


Clearly Ep{Hq) is finite. Moreover for any n > 1 

\dHi^ I < ((A„ © B{0, R)) \ (A„ © B{0, R))) < AdpR{2n + 2R) 


d-l 


which implies [MeanEnergy]. Similar computations show that assumption [Boundary] holds 
as well. Finally, Assumption [Regularity] is a consequence of the following uniform continuity 
of Hq coming from a simple geometric analysis: for any 9 = {z,R,f3), 9' = {z',R',(i') in 0 (we 
assume R< R') 


\H^o - <1 < 1^0 




-R, 


C'.A/3)| 




\H, 


(z\R',P) 


-R, 




<\z- z’\Na, + 2dm2R'r\R - R'\Namo,r') + I/? - /3'|. 


5.7 Proof of Proposition [2] 

Since the Quermass-interaction model has an exponential form, it is sufficient to check the assump¬ 
tions of Corollary[TJ The existence of the model for any z > 0 and any /3 G is proved in Q while 
stability of the interaction follows from [1^. The assumptions [MeanEnergy] and [Boundary] 
can be proved as in Section [foU see (l64ll . thanks to the additivity of the Minkowski functionals. 
The assumption [VarPrin] follows from y since the interaction is stable and finite range. 
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